Excitation energies, polarizabilities, multipole transition rates, and lifetimes in Th IV 
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Excitation energies of the ns 1 / 2 (n= 7-10), npj (n = 7-9), ndj (n= 6-8), nfj (n= 5-7), and ngj 
(n— 5-6) states in Th IV are evaluated. First-, second-, third-, and all-order Coulomb energies and 
first- and second-order Coulomb-Breit energies are calculated. Reduced matrix elements, oscillator 
strengths, transition rates, and lifetimes are determined for the 96 possible nlj —n'lji electric-dipole 
transitions. Multipole matrix elements (7si/2 — 6dj,7s 1 / 2 — 5/j , and 5/5/2 — 5/7/2) are evaluated 
to obtain the lifetimes of the 5/7/2 and 7si/2 states. Matrix elements are calculated using both 
relativistic many-body perturbation theory, complete through third order, and a relativistic all- 
order method restricted to single and double (SD) excitations. Scalar and tensor polarizabilities for 
the 5/ 5 /2 ground state in Th + are calculated using relativistic third-order and all-order methods. 
These calculations provide a theoretical benchmark for comparison with experiment and theory. 

PACS numbers: 31.15.Ar, 31.15.Md, 32.10.Fn, 32.70.Cs 



I. INTRODUCTION 



A detailed investigation of radiative parameters for 
electric dipole (El) transitions in Fr-likc ions with Z 
= 89-92 was presented recently by Biemont et al. 0, 
where relativistic Hartree-Fock and Dirac-Fock atomic 
structure codes were used to calculate transition rates 
and oscillator strengths for a limited number of tran- 
sitions using experimental energies given by Blaise and 
Wyart [2|. In the compilation [2|, experimental energies 
are given for 56 levels of neutral Fr, 24 levels of Fr-like 
Th, and seven levels of Fr-like Ac and U. Experimental 
energies for 13 levels of Fr-like Ra were reported in the 
NIST compilation 0. 

Lifetime measurements for the 7pj, 6dj, 9sj/2, and 
8si/2levels of neutral francium were presented in Refs. |4j, 
01 El lI HI < In those papers, experimental measurements 
were compared with ab-initio calculations performed by 
Johnson et al. @, by Dzuba et al. [13,113, by Safronova 
et al. [13, and by Safr onova and Johnson [13j|. Third- 
order many-body perturbation theory was used in Ref. 
to obtain El transition amplitudes for neutral alkali- 
metal atoms. The correlation potential method and the 
Feynman diagram technique was used in Refs. [IfllTlT ] to 
calculate El dipole matrix elements in neutral francium 
and in Fr-likc radium. Calculations of atomic proper- 
ties of alkali-metal atoms in Refs. fl2l Hal were based 
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on the relativistic single-double (SD) approximation in 
which single and double excitations of Dirac-Fock wave 
functions were included in all orders in perturbation the- 
ory. 

In the present paper, relativistic many-body pertur- 
bation theory (RMBPT) is used to determine energies, 
matrix elements, oscillator strengths, and transition rates 
for multipole transitions in Fr-like thorium. These cal- 
culations start from a radonlike closed-shell Dirac-Fock 
(DF) potential. It should be noted that Th IV is the first 
ion in francium isoelectronic sequence with a [Rn] 5/5/2 
ground state instead of a [Rn]7si « ground state, as for 
Fr I, Ra II, and As III. Correlation corrections become 
very important for such systems as was recently demon- 
strated in Savukov et al. |14| . where it was shown that 
the ratio of the second-order to lowest-order removal en- 
ergy for the [Xe] 4/5/2 ground state in Ce IV and Pr V is 
18 % and 11 %, respectively. 

We calculate excitation energies of nsi/2 (n= 7-10), 
npj (n = 7-9), ndj (n= 6-8), nfj (n— 5-7), and ngj (n= 
5-6) states in Fr-like thorium. Reduced matrix elements, 
oscillator strengths, transition rates, and lifetimes are de- 
termined for the 96 possible nlj — n'l'j, electric-dipole 
transitions. Multipole matrix elements (7s!/ 2 — 6dj, 
7si/2 — 5/j, and 5/ 5 / 2 — 5/7/2) are evaluated to ob- 
tain the lifetimes of 5/7/2 and 7si/2 states. Scalar and 
tensor polarizabilities of the 5/5/2 ground state of Th 3+ 
are also calculated. Matrix elements are calculated using 
both relativistic many-body perturbation theory, com- 
plete through third order, and the relativistic all-order 
method restricted to single and double (SD) excitations. 
Such calculations permit one to investigate the conver- 
gence of perturbation theory and estimate the error in 
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theoretical data. 



II. THIRD-ORDER AND ALL-ORDER RMBPT 
CALCULATIONS OF ENERGIES 

As mentioned in the introduction, we carry out all of 
the calculations in this work using two methods, third- 
order MBPT, described in 0, and the relativistic all- 
order SD method, described in 0, and references 
therein. The SD method includes correlation corrections 
in a more complete way and is expected to yield more 
accurate results, especially when correlation corrections 
are significant. While the SD method includes fourth- 
and higher-order terms, it omits some third-order terms. 
These omitted terms are identified and added to our SD 
data (see 16] for details). 

We use the B-spline method ^7] to generate a com- 
plete set of basis DF orbitals for use in the evaluation 
of RMBPT expressions. For Th IV, we use 50 splines 
of order k = 8 for each angular momentum. The basis 
orbitals are constrained to a spherical cavity of radius R 
— 45 a.u.. The cavity radius is chosen large enough to 
accommodate all nlj orbitals considered here and small 
enough that 50 splines can approximate inner-shell DF 
wave functions with good precision. We use only 40 out of 
50 basis orbitals for each partial wave in our third-order 
and all-order energy calculations since contributions from 
higher-energy orbitals are negligible. 

Results of our energy calculations for the 28 states 
of Th IV are summarized in Table d Columns 2-7 of 
Table d give the lowest-order DF energies E^-°\ second- 
and third-order Coulomb correlation energies, and 
first-order Breit contribution B^\ second-order 
Coulomb-Breit corrections B^ 2 ', and the Lamb shift £ls- 
The sum of these six contributions is our final third-order 
RMBPT result E ( t H listed in the eighth column of Tabled 
The all-order SD energies are listed in the column E ST> , 
and that part of the third-order energy omitted in the 

(3) 

SD calculation is given in column -E^xtra- We note that 
E SD includes E^ completely. We take the sum of the 
six terms E<°\ E SD , E^ tIa , B {1 \ B^ 2 \ and E LS to be 
our final all-order result Ef^, listed in the eleventh col- 
umn of Table d Experimental energies from Blaise and 
Wyart @ are given in the column labeled E cxpt . Dif- 
ferences between third-order and experimental energies 
(5_E(3) = fiW _ E C x P t, and between SD and experimental 
energies 5E SU — Ef^ — E cxpt , are given in the last two 
columns of Table d respectively. 

As expected, the largest correlation contribution to the 
valence energy comes from the second-order term E^ 2 \ 
This term is relatively simple to calculate; thus, we cal- 
culate E^ with better numerical accuracy than E^ 
and E SD . The second-order energy E^ includes par- 
tial waves up to Z max = 8 and is extrapolated to account 
for contributions from higher partial waves (see, for ex- 
ample, [l8L IT^|). As an example of the convergence of 
E^ with the number of partial waves I, consider the 



5/5/2 state in Th IV. Calculations of E^ with Z max = 

6 and 8 yield £ (2) (5/ 5 / 2 ) = -30810 and -31771 cm" 1 , 
respectively. Extrapolation of these calculations yields 
—32100 and —32154 cm -1 , respectively. Therefore, we 
estimate the numerical uncertainty of E^ (5/5/2) to be 
approximately 54 cm . It should be noted that this is 
the largest contribution from the higher partial waves, 
since we obtain a numerical uncertainty of 26 cm -1 for 
E^^dj) and the numerical uncertainty of 1 cm -1 for 
£ (2) (7sj). Similar convergence patterns are found for all 
other states considered. 

We use Z max = 6 in our all-order calculations owing to 
the numerical complexity of the E su calculation. As we 
noted above, the second-order E^ is contained in the 
_B SD value. Therefore, we use our high-precision calcu- 
lation of E^ described above to account for the contri- 
butions of the higher partial waves. We simply replace 

[Zmax = 6] value with the final high-precision second- 
ly) 

order value -Ei The same number of partial waves, 
'max = 6, is used in the third-order calculation. Since 
the asymptotic Z-dependence of the second- and third- 
order energies are similar (both fall off as l~ 4 ), we use 
the second-order remainder as a guide to estimate the 
numerical errors in the third-order contribution. The 
contribution -E^xtra given in Table d accounts for that 
part of the third-order RMBPT correction not included 

(3) 

in the SD energy. The values of E^xtra are smaller than 
the values of E^> by approximately a factor of 3. 

The first-order Breit energies (column #W of Tabled) 
include retardation, whereas the second-order Coulomb- 
Breit energies (column B<-V of Tabled) are evaluated us- 

(3) 

ing the unretarded Breit operator. The total E tot in 
Table ID is the sum of six terms, E^°\ £ (2) , £ (3) , 
B^ 2 \ and -Els- We find that the correlation corrections 
to energies are large, especially for the 5/j states. For 
example, E^ is about 15% of E^ and E^ is about 
36% of E^ for the 5/j states. Despite the evident slow 
convergence of the perturbation theory expansion, the 
5/j energy from the third-order RMBPT calculation is 
within 0.9% of the measured energy. It should be noted 
that correlation corrections are much smaller for all other 
states; the ratios of E^ and E^ are equal to 6%, 5%, 
and 2% for the 6dj, 7sx/2, and 10sx/2 states, respectively. 
An important consequence of the large size of correlation 
corrections for 5/j states is a different ordering of un- 
corrected and correlated energies. As can be seen from 
Tabled —E^ values for 6dj states are larger than —E^ 
values for 5/j states; however, — (£^°) + E^) values for 
6dj states are smaller than — (£^°) + E^) values for 5/j 
states; thus, although DF calculations predict the ground 
state of Th IV to be 6d 3 / 2 , correlated calculations cor- 
rectly predict the ground state to be 5/ 5 / 2 . 

The quantity E s £ t in Tabled is the sum of six terms; 
j 5 (o) j ^sd^ s (3) raj B (2) j and e ls)_ The column 

labeled 5E SU in Tabled gives differences between our ab 
initio results and the experimental values 2\. The SD 
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TABLE I: Zeroth- (DF), second-, and third-order Coulomb correlation energies E^ n \ single-double Coulomb energies E SD , 
.Egxtrai first-order Breit and second-order Coulomb-Breit corrections to the energies of Th IV. The total energies {E^ t = 
E (0) +E (2) +E (3) +B (i) +B (2) + £ LS) E fD = E (o) + E SD +E^+B m +_b(2) + Els ) for Th IV are compared with experimental 
energies £ exp t 0, SE = E to t - E c * p t. Units: cm -1 . 



nlj 


S (0) 


E (2) 


E (3) 


5(1) 


5(2) 


E LS 


g(3) 
a tot 


E su 


E (3) 

extra 


E tot 


-Ecxpt 


8E^ 


SE SD 


5/5/2 


-206606 


-32100 


11739 


704 


-2747 





-229010 


-26327 


4672 


-230304 


-231065 


2055 


761 


5/7/2 


-203182 


-30549 


10954 


521 


-2616 





-224872 


-25252 


4361 


-226168 


-226740 


1868 


571 


6d3/2 


-211799 


-13258 


4129 


438 


-880 





-221370 


-11422 


1663 


-222000 


-221872 


502 


-128 


6d5/2 


-207574 


-11608 


3300 


326 


-807 





-216364 


-10208 


1337 


-216927 


-216579 


214 


-348 


7S1/2 


-200273 


-11204 


4402 


325 


-458 


89 


-207119 


-9455 


1697 


-208075 


-207934 


815 


-140 


7P1/2 


-165095 


-7991 


2782 


298 


-272 





-170278 


-7147 


1125 


-171091 


-170826 


548 


-265 


7f>3/2 


-153572 


-6213 


2124 


184 


-226 


1 


-157703 


-5619 


861 


-158372 


-158009 


306 


-363 


8S1/2 


-109201 


-3840 


1462 


122 


-163 


25 


-111594 


-3255 


575 


-111896 


-111443 


-151 


-453 


7d 3 /2 


-108639 


-3815 


1110 


113 


-199 





-111429 


-3559 


496 


-111787 


-111380 


-49 


-406 


7d 5/2 


-107032 


-3592 


1011 


86 


-188 





-109715 


-3401 


452 


-110083 


-109638 


-77 


-445 


6/5/2 


-99921 


-5078 


1615 


77 


-270 





-103577 


-4763 


710 


-104167 


-103796 


219 


-371 


6/7/2 


-99481 


-4998 


1590 


61 


-276 





-103103 


-4644 


696 


-103644 


-103250 


146 


-394 


8P1/2 


-94597 


-3199 


1112 


129 


-116 





-96670 


-2861 


460 


-96984 


-96549 


-122 


-436 


8P3/2 


-89595 


-2628 


898 


83 


-101 





-91342 


-2397 


371 


-91638 


-91194 


-148 


-444 


5fl7/2 


-70583 


-1506 


421 





-4 





-71672 


-1592 


207 


-71971 


-71694 


22 


-296 


5ff9/2 


-70602 


-1491 


422 





-4 





-71675 


-1548 


205 


-71949 


-71675 





-255 


9si /2 


-69455 


-1849 


684 


61 


-80 


11 


-70628 


-1620 


275 


-70810 


-70337 


-291 


-473 


8d 3 / 2 


-68677 


-1846 


522 


54 


-93 





-70040 


-1821 


244 


-70292 


-69111 


-929 


-1181 


8d 5 / 2 


-67864 


-1767 


488 


42 


-89 





-69190 


-1780 


228 


-69463 


-68537 


-653 


-926 


7/5/2 


-64674 


-2576 


798 


40 


-137 





-66549 


-2504 


362 


-66912 


-66510 


-39 


-402 


7/7/2 


-64454 


-2517 


782 


31 


-139 





-66296 


-2409 


352 


-66619 


-66006 


-290 


-613 


9P1/2 


-61858 


-1663 


566 


68 


-61 





-62948 


-1548 


241 


-63158 








9P3/2 


-59200 


-1400 


468 


45 


-54 





-60141 


-1352 


199 


-60362 








6fl7/2 


-49121 


-1042 


281 





-4 





-49884 


-1242 


145 


-50221 


-49391 


-493 


-831 


6fl9/2 


-49138 


-1022 


285 





-4 





-49879 


-1166 


142 


-50166 


-49390 


-489 


-775 


10ai/a 


-48177 


-1046 


376 


35 


-46 


5 


-48853 


-990 


154 


-49018 


-48624 


-229 


-394 



results agree better with the experimental values than 
the third-order RMBPT results for low-lying states where 
the correlation correction is larger. Comparison of the 
results given in two last columns of Table [I] shows that 
the ratio of SE^ and SE SD is about 3 for the 5fj states. 
As expected, including correlation to all orders led to 
significant improvement of the results. Better agreement 
of the all-order values with experiment demonstrates the 
importance of the higher order correlation contributions. 



III. ELECTRIC-DIPOLE MATRIX ELEMENTS, 
OSCILLATOR STRENGTHS, TRANSITION 
RATES, AND LIFETIMES IN TH IV 

A. Electric-dipole matrix elements 

The calculation of the transition matrix elements pro- 
vides another test of the quality of atomic-structure cal- 
culations and another measure of the size of correlation 
corrections. Reduced electric-dipole matrix elements be- 
tween low-lying states of Th IV calculated in third order 
RMBPT and in the SD approximation are presented in 
Table |nl 

Third-order matrix elements Z^ DF+2+3 ^ include DF 



contributions together with second-order Z^ 2 > and third- 
order Z^ 1 correlation corrections. Second- and third- 
order random-phase-approximation (RPA) terms are it- 
erated to all orders in the present calculation. Third- 
order corrections include Brueckner orbital (BO), struc- 
tural radiation Z^ SRS> , and normalization ^( NORM ) cor- 
rections, in addition to the third-order RPA terms, see 
0. The terms Z ( RPA ) and Z ( BO ) give the largest con- 
tributions to the total. The sum of terms Z ( RPA ) and 
Z ( B °) is about 15-25% of Z( DF > and has a different 
sign. Structural radiation, and normalization corrections 
are small. We find correlation corrections Z^ 2+3 ^ to be 
very large, 10-25%, for many cases. All results given 
in Table ^ arc obtained using length-form matrix ele- 
ments. Length-form and velocity-form matrix elements 
differ typically by 5-20 % for DF matrix elements and 
2-5 % for the second-order matrix elements in these cal- 
culations. 

Electric-dipole matrix elements evaluated in the all- 
order SD approximation are given in columns labeled 
Z( SD ) of Table |TT| The SD matrix elements Z^ in- 
clude Z^ completely, along with important fourth- and 
higher-order corrections. The fourth-order corrections 
omitted from the SD matrix elements were discussed re- 
cently by Derevianko and Emmons poj . The SD matrix 
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TABLE II: Reduced electric-dipole matrix elements calculated to first, second, third, and all orders of RMBPT in Th IV. 



Transition 


^T^ur ) 






Transition 


r-7(TtW\ 

^\ uc > 






Fifl-T In 

u tJ7/2 


Ft ft: In 


1.1236 


0.6400 


0.7118 


(\d C In 

uu, 5/2 


VJ7/2 


3.3539 


3.1272 


3.0192 


^97/2 


Fi 1*7 In 
J J 7/2 


0.2298 


0.1344 


0.1511 


6l^r In 
WUj o / 2 


' «/ 5/2 


0.3586 


0.2561 


0.2398 


^97/2 


B T r In 

5/2 


9.9381 


8.8278 


8.8815 


6rtr In 
wtA 5 / 2 


7/7/2 


1.5854 


1.1995 


1.1083 


Bflfl In 


6 f*f in 
^ J 7/2 


1.9244 


1.7141 


1.7252 


V)fJ,r In 


1 i J Aj 2 


3.1975 


2.7006 


2.7549 


^97/2 


' J5/2 


8.5105 


8.6350 


8.4527 


BWr in 

uu, 5/2 


°f 3/2 


0.6529 


0.4291 


0.4288 


tiO-v In 

a y7/ 2 


7/7/2 


1.6217 


1.6446 


1.6082 


6rtr In 


9eJq /o 
ty /-'o / 2 


0.3496 


0.1743 


0.1438 


in 


~V 5/2 


0.8677 


0.4157 


0.4134 


1 ^5 / 2 


fr /o 

U J 5/2 


0.0048 


0.0360 


0.0521 


(\fl-7 In 


PS f-7 In 

°J 7/2 


0.1757 


0.0766 


0.0918 


' ^5/2 


Ft Tt In 

°J 7/2 


0.0703 


0.2298 


0.2961 


u i/7/2 


fS fr in 


1.3517 


1.2696 


1.1321 


' ^5/2 


(\fr In 

u «/5/2 


2.0852 


1.8440 


1.8605 


^^7/2 


fi fv /o 

u J 7/2 


0.2464 


0.2319 


0.2042 


' ^5 / 2 


f 7 
u .y 7/2 


9.3944 


8.3249 


8.4004 


yj yj/2 


1 J 5/2 


14.6139 


13.4006 


13.5178 


7dr in 
' u '5/2 


' J 5/2 


1.0359 


1.1377 


1.1176 


u y7/ 2 


7/7/2 


2.8402 


2.6137 


2.6402 


7d<r in 
' ^5/ 2 


7/7/2 


4.5090 


4.9585 


4.8615 












7fl,r in 
' ^0 / 2 


' FS/2 


5.9481 


5.3934 


5.4192 




7/2 


1.3635 


0.8187 


0.9054 


7^5/2 


°f 3/2 


6.8642 


6.5224 


6.5180 


5<?9/2 


6/7/2 


11.3858 


10.1422 


10.2068 


7c?5/2 


9p3/2 


0.8580 


0.7746 


0.7645 


^^9/2 


7/7/2 


9.5778 


9.7291 


9.5435 












u y9/2 


Ft It In 

°J 7/2 


1.0425 


0.5803 


0.5682 


Rdr in 


5fr /O 


0.0062 


0.0307 


0.0461 


6^9/2 


6/7/2 


1.4474 


1.3738 


1.2252 


8r/r ,n 


5/7/2 


0.0082 


0.2138 


0.2520 


6^9/2 


7f^ /« 
' 77/2 


16.8152 


15.4636 


15.5953 


Rrlc In 


f c /o 
u 75/2 


0.3783 


0.3187 


0.3357 












OLt 5/2 


u l / 7/2 


1.7787 


1.5158 


1.6020 


(Sf/o /n 


dj5/2 


2.4281 


1.3367 


1.5295 


8c?5/2 


1 J 5/2 


3.5961 


3.3715 


3.3634 


(\d,o in 


U J 5/2 


2.6761 


2.4407 


2.3610 


8c?5/2 


7/7/2 


16.1550 


15.1766 


15.1464 


(\fir> In 
ULt 3/2 


' 75/2 


1.2888 


0.9082 


0.8469 


OLt 5/2 


1 f 3/2 


1.1741 


0.8883 


0.8651 


(\fir> In 


' Fl/2 


2.5465 


2.0723 


2.1220 


OLt 5/2 


°f 3/2 


9.3341 


8.9653 


8.9992 


6^3/2 


' K3/2 


0.9963 


0.8270 


0.8488 


8(Ve /o 

OU-5/2 


&PS/2 


11.4942 


11.1245 


11.0183 


6(^3/2 


8pi/2 


0.4074 


0.1809 


0.1906 












(\(ir> In 


3/2 


0.2173 


0.1410 


0.1422 


7si /-> 




2.8994 


2.3669 


2.4196 






0.2148 


0.0314 


0.0221 


7Si/2 


7P3/2 


3.9933 


3.2930 


3.3677 


6d 3 / 2 


9P3/2 


0.1172 


0.0568 


0.0475 


7S1/2 


8P1/2 


0.0565 


0.2438 


0.2346 








7si /o 


°f 3/2 


0.3273 


0.0639 


0.0679 


' ti 3/2 


3/5/2 


0.0654 


0.2134 


0.2587 


7cii In 




0.0550 


0.1958 


0.1849 


' ti 3/2 


6/5/2 


7.8264 


6.9261 


6.9826 


' *l/2 


y P3/2 


0.1326 


0.0706 


0.0674 


' ti 3/2 


7f c /o 
' J5/2 


3.5114 


3.8940 


3.8365 












7 do In 


7t>i 

l Pl/2 


3.8261 


3.4234 


3.4490 


oa l/2 


1/2 


1.5874 


1.5601 


1.5492 


Id-i in 
1 ^3 / 2 


7 Do in 


2.0308 


1.8374 


1.8444 


8.S1 /o 
1/2 


7wo /o 


3.0768 


3.0010 


2.9756 


7 do in 
1 ^3/2 


ftiT)-i In 


5.4788 


5.1788 


5.1791 


8^1 


°f 1/2 


5.0325 


4.6912 


4.7280 


1 ^3/2 


HT)o /o 

°K3/2 


2.1716 


2.0607 


2.0630 


8^1 

°^l/2 


°r 3/2 


6.7737 


6.3299 


6.3880 


1 Ll vJ / 2 


1/2 


0.4173 


0.3208 


0.3214 


8,Si /o 

1/2 


1/2 


0.0758 


0.1951 


0.2099 


» ^3/2 


9p3/2 


0.3057 


0.2822 


0.2790 


8^1/2 


9p 3 /2 


0.5353 


0.4009 


0.3719 


o&3/2 


3J5/2 


0.0029 


0.1874 


0.2098 




' Pl/2 


0.4722 


0.4773 


0.4657 




6/5/2 


1.6194 


1.3910 


1.4493 


9S1/2 


7P3/2 


0.7567 


0.7254 


0.7123 


8d 3/2 


7/5/2 


13.4659 


12.6494 


12.6151 


9S1/2 


8P1/2 


2.8858 


2.8094 


2.8070 


8^3/2 


7P1/2 


0.9598 


0.7550 


0.7418 


9S1/2 


8J53/2 


5.4176 


5.3126 


5.2854 


8^3/2 


7p 3 /2 


0.3604 


0.2581 


0.2523 


9Si/2 


9P1/2 


7.6975 


7.4191 


7.4316 


8g! 3 /2 


8P1/2 


5.9596 


5.6802 


5.7284 


9si/ 2 


9p 3 / 2 


10.2453 


9.8671 


9.8961 


8g! 3 /2 


8P3/2 


3.2096 


3.0818 


3.0900 












8g! 3 /2 


9P1/2 


9.1563 


8.8665 


8.8044 


lOsi/2 


7P1/2 


0.2605 


0.2651 


0.2489 


8rf 3 /2 


9P3/2 


3.6602 


3.5363 


3.5120 


IOS1/2 


7P3/2 


0.4011 


0.3784 


0.3596 












IOS1/2 


8P1/2 


0.7932 


0.7936 


0.7873 


6d 5 / 2 


5/5/2 


0.6391 


0.3624 


0.4116 


10si/2 


8P3/2 


1.2032 


1.1808 


1.1676 


6o! 5 / 2 


5/7/2 


2.9557 


1.7032 


1.9190 


IOS1/2 


9P1/2 


4.5019 


4.3611 


4.3523 


6^5/2 


6/5/2 


0.7669 


0.7085 


0.6847 


10Sl/2 


9P3/2 


8.3091 


8.1560 


8.0743 
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TABLE III: Comparison of length [L] and velocity [V] results 
for reduced electric-dipole matrix elements in lowest and third 
orders of perturbation theory in Th IV. 
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1 


7779 


8d 5 /2 


5/7/2 


-0 


0082 


-0 
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3256 


7si /2 
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8994 


2 


6904 


2 


3943 


2 


3936 
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7P3/2 


-3 


9933 


-3 


7030 


-3 


3434 


-3 


3432 


7Si/ 2 


8P1/2 


-0 


0565 


-0 


0910 


-0 


2334 


-0 


2335 


8*1/2 


9P1/2 


-0 


0758 


-0 
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-0 


1931 


-0 


1932 


8S1/2 
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4031 
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-0 


4722 


-0 


4420 


-0 


4737 


-0 


4736 
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6966 





7181 





7180 


9*1/2 


8P1/2 


-2 


8858 


-2 


8249 


-2 


8139 


-2 


8137 



elements ij( SD ) are smaller than Z( DF+2 ) but larger than 
z (df+2+3) for all of the trans itions listed in Table |TT| 



B. Form-independent third-order transition 
amplitudes 

We calculate electric-dipole reduced matrix elements 
using the form-independent third-order perturbation the- 
ory developed by Savukov and Johnson in Ref . [2l[ . The 
precision of this method has been demonstrated pre- 
viously for alkali-metal atoms. In this method, form- 
dependent "bare" amplitudes are replaced with form- 
independent random-phase approximation ( "dressed" ) 
amplitudes to obtain form-independent third-order am- 
plitudes. As in the case of the third-order energy calcu- 
lation, a limited number of partial waves with Z max < 7 
is included, giving rise to some loss of gauge invariance. 
Comparison of length- and velocity-form matrix elements 
serves as a measure of the numerical accuracy of the re- 
sulting calculations. 

Length- and velocity-form matrix elements from DF, 
second-order, and third-order calculations are given in 
Table IIIII for the limited number transitions in Th IV. 
Following the procedure discussed in Ref. j^, the DF 
and RPA matrix elements in the table were obtained by 
dividing the corresponding amplitude by the lowest-order 
transition energies while the third-order matrix elements 
were obtained by dividing the third-order amplitude by 
the second-order transition energies. Values of Z^ F > dif- 



TABLE IV: Wavelengths A (A), weighted transition rates 
gA (s _1 ), and oscillator strengths gf in Th IV. The SD data 
(g,4 (SD) and ff/ (SD) ) are compared with theoretical (gA (iIFR) 
and <?/^ HFR ') values given in Ref. Numbers in brackets 
represent powers of 10. 
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+0] 


1.1 [+0] 


5/5/2 


637/2 
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TABLE V: Lifetimes r in ns for the nl levels in Fr-like Th. 


Level 




r (SD) 




Level 






r (SD) 
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7P1/2 






1.099 


6d 5 /2 




676. 




7P3/2 






0.632 


7d 3 / 2 




0.667 




8P1/2 






3.194 


7d 5 /2 




0.854 




8P3/2 






1.871 


8e?3/2 




1.176 




9P1/2 






5.893 


8(^5/2 




1.600 




9P3/2 






4.933 


537/2 




0.815 




6/5/2 






0.300 


5ff9/2 




0.780 




6/7/2 






0.297 


637/2 




1.768 




7/5/2 






0.684 


639/2 




1.567 




7/7/2 






0.639 


8S1/2 
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fer in L and V forms by 2-15% for the p-s transitions. 
Huge L — V differences in the Z( DF ) for d-f transitions 
can be seen in Table IIIII Third-order calculations es- 
sentially remove such differences; the residual differences 
(0.002%-0.2%) being explained by the limited number of 
partial waves used in the evaluation of third-order matrix 
elements. 
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C. Oscillator strengths, transition rates and 
lifetimes 



We calculate oscillator strengths and transition prob- 
abilities for 96 possible nlj — n'l'j, electric-dipole tran- 
sitions including the ns\/2 (n= 7-10), npj (n — 7-9), 
ndj (n= 6-8), nfj (n— 5-7), and ngj (n= 5-6) states in 
Fr-like thorium. Our results are presented in Table IIVI 
and Table E| Wavelengths A (A), weighted transition 
rates gA (s^ 1 ), and oscillator strengths gf in Th IV are 
given in Table llVl Our SD data, 5 ^ SD ) and 5/ (SD) , 
are compared with theoretical calculations, gA^ RFR ^ and 
g/ (HFR) , from Ref. Q. It should be noted that experi- 
mental energies are used to calculate gA^ SD ^ and gf^ SD ^ 
as well as gA^ HFR ^ and gf^ HFR ^- Therefore, we really 
compare the dipole matrix elements (see Table ITTf) . The 
SD and HFR results for s — p and p — d transitions dis- 
agree by 6-25%, except for the 7pi/ 2 — 8s 1 / 2 transition 
with 60% disagreement. There are also substantial dis- 
agreements (factors of 2-5) between SD and HFR results 
for the / — g and / — d transitions. Correlation correc- 
tions are very important for those transitions as discussed 
above (see Table ITTjl . The RPA and BO contributions 
have the same sign opposite to the DF contributions and 
the total values are half of the DF values. We see from 
Table llVl that for the f — g and / — d transitions the val- 
ues of gA^ HFR ^ and <?/( HFR ) are larger by a factor of 2-5 
than the gA^ SD ^ and gf^ SD ^ values, respectively. On the 
basis of these comparisons, it appears that correlation 
corrections were not included in Ref. [| for transitions 
involving the 5/j states. Our conclusion is confirmed 
by comparison of gA^ RFR ^ and g/ < - HFR - ) with our gA^° F ^ 
and gf( DF ^ results (compare the Z DF and Z SD columns 
in Table [HJ). The disagreement between HFR and DF 
values for transitions rates and oscillator strengths is sig- 
nificantly smaller (about 10-20 %) than the disagreement 
between HFR and SD (by a factor of 2-5). 



We calculate lifetimes of ns\i 2 (n= 8-10), npj (n = 7- 
9), ndj (n— 6-8), nfj (n— 6-7), and ngj (n= 5-6) states 
in Fr-like thorium using the SD results for dipole matrix 
elements and experimental energies Q- We list lifetimes 
T (SD) m Xable \V\ Unfortunately, there are no experi- 
mental measurements to compare with our results; how- 
ever, we hope that our calculations provide a theoreti- 
cal benchmark and lifetime measurements will be carried 
out. 



Lifetimes for two excited levels, 7si/2 and 5/7/2 were 
not included in Table since there are no electric- 
dipole transitions from these levels. Contributions of the 
electric- and magnetic-multipole transitions to the life- 
time of the 7si/2 and 5/7/2 levels are considered below. 



IV. MULTIPOLE MATRIX ELEMENTS, 
TRANSITION RATES, AND LIFETIMES IN 
TH IV 



Reduced matrix elements of the electric-quadrupole 
(E2), electric-octupole (E3), and magnetic-multipole 
(Ml, M2, and M3) operators in lowest, second, third, and 
all orders of perturbation theory are given in Table I VII 
for Th IV. Detailed description of the calculations of the 
multipole matrix elements in lowest and second orders 
of perturbation theory were given in Refs. [22I I23I |24| . 
Third-order and all-order calculations are done in the 
same way as the calculations of the El matrix elements. 
In Table IVT1 we present E2, E3, Ml, M2, and M3 matrix 
elements in the Z^ F \ Z( DF+2 \ z( DF+2 + 3 ), and 
approximations for the 5/ 5 / 2 — 5/7/2, 5/j — 7s X j 2 , and 
6dj — 7s!/ 2 transitions in Th IV. 

The second-order contribution is about 1-3% for all 
transitions involving the 7s!/ 2 states, but it is different 
for the 5/5/2 — 5/7/2 transition. It is very small (0.1%) for 
the Ml matrix elements and rather large (20%) for the 
E2 matrix elements. The large difference between Z^° F ' 
and Z (DF+2+ V or Z( SD > for E2 and E3 operators could be 
explained by the large size of the Brueckner orbital (BO) 
correction; the ratios of Z^ B °^ and Z^° F ^ are equal to 0.06 
and 0.3 for the 6dj — 7s x / 2 and 5fj — 7s!/ 2 transitions, 
respectively. 

Wavelengths and transition rates A( SD ) for the electric 
multipole (E2 and E3) and magnetic-multipole (Ml, M2, 
and M3) transitions in Th IV calculated in the SD ap- 
proximation are presented in Table IV1TI The largest con- 
tribution to the lifetime of the 5/7/2 state comes from the 
Ml transition. The largest contribution to the lifetime of 
the 7si/2 state comes from the E2 transitions. Our SD 
result for Ml matrix elements are in perfect agreement 
(0.5%) with HFR results obtained by by Biemont et al. 
111 Ref. 0. The disagreement is much larger between 
HFR and SD results for E2 transitions; 6% and 18% for 
the 6g? 5 /2 — 7s!/ 2 and 6d 5 / 2 — 7s l / 2 matrix elements, re- 
spectively This is expected because the correlation cor- 
rection to the Ml matrix clement is very small, but the 
correlation correction to E2 matrix element is large, as 
discussed above. Since the SD all-order method includes 
the correlation corrections in a rather complete way, we 
expect to see disagreements with HFR calculations in the 
cases where correlation corrections are significant. 

Finally, we find that the lifetime of the 5/7/2 state is 
1.07 s and the lifetime of the 7si/2 state is 0.590 s. An 
estimate of the 7s!/ 2 state lifetime (about 1 s) was given 
by Peik and Tamm in Ref. |25| . but no measurements 
have as yet been performed. 
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TABLE VI: Reduced matrix elements of the electric multipole (E2, E3) and magnetic-multipole (Ml, M2, and M3) operators 
in first, second, third, and all orders of perturbation theory in Th IV. 



Transition £< DF+2 > ^(df+2+3) 

Ml 5/5/2 5/7/2 1-8506 1.8525 1.8390 1.8514 

E2 5/5/2 5/ 7/2 1.5669 1.2339 0.9724 1.0834 

E2 6d 3/2 7s 1/2 7.7806 7.5518 6.9232 7.0631 

E2 6d 5/2 7s 1/2 10.0084 9.7894 8.9992 9.1526 

E3 5/ 5/2 7s 1/2 12.9552 13.1718 9.3847 10.6349 

E3 5/ 7/2 7s 1/2 15.7681 16.0030 11.5616 13.0648 

M3 6d 5/2 7si /2 59.2100 60.9891 59.7136 57.3087 



TABLE VII: Wavelengths A (A) and transition rates A (SD) 
(s _1 ) of the electric multipole (E2, E3) and magnetic- 
multipole Ml, M2, and M3) transitions in Th IV calculated in 
the SD approximation. Numbers in brackets represent powers 
of 10. 





Transition 


A 


^(SD) 


Ml 


5/5/2 


5/7/2 


23119.6 


9.352[-l] 


E2 


5/5/2 


5/7/2 


23119.6 


2.487[-5] 


E2 


6d 3 / 2 


7si /2 


7174.88 


1.469[-0] 


E2 


6d 5 /2 


7si /2 


11568.2 


2.264[-l] 


Ml 


6^3/2 


7si /2 


7174.88 


5.429[-5] 


E3 


5/5/2 


7si/ 2 


4323.25 


6.299[-7] 


E3 


5/7/2 


7si /2 


5317.62 


2.232[-7] 


M2 


5/5/2 


7si /2 


4323.25 


5.154[-9] 


M3 


6d 5 /2 


7si/ 2 


11568.2 


1.863 [-8] 



V. SCALAR AND TENSOR 
POLARIZABILITIES IN THE 5/ 5/2 GROUND 
STATE OF TH 3 + 

We calculate the tensor polarizability q 2 of Th 3+ in a 
state v using a sum-over-states approach |26| 



a 2 (v) = ^I(v 



(1) 



Here 



I(v — n) = A(vn)- 



7 2 

vn 



(2) 



where 



A(nv) = 



5j v (2jv - 1) 



y 6(j v + l)(2j„ + l)(2j v + 3) 



jv 1 jn 

1 j v 2 



is a reduced electric-dipole matrix element. 
The calculation of the 0:2(5/5/2) is divided into three 



parts: 

a2 lain (5/5/2) 

SaT Ie (5/5/2) 

Of' (5/5/2) 



^2 7 ( 5 /5/2 - ndj) + ^2 1(5 f 5/2 - ng 7/2 ); 

n— 6 
5 

^ J ( 5 /5/2 - ndj); 

n— 3 
50 

J ( 5 /5/2 - ndj) + ^2 1(5/5/2 - ng 7/2 () 



n— 5 



50 



n=9 



n=7 



We present the details our calculations of tensor po- 
larizabilities a 2 for the ground state 5/5/2 in Table I*VTTT1 
We use experimental energies from 2j. Electric-dipole 
matrix elements evaluated in the the third-order and all- 
order SD approximations are given in columns labeled 
^(DF+2+3) anc j ^(sd)^ correS p 0nc ii n g contributions 

to the tensor polarizability are given in columns labeled 

(DF4 
*2 



^ and a 2 SD ^ 



o 2 ail and Sa c 2 ° 



o 2 - ana a 2 ' . The contributions 

are found to be very small and are calculated in the DF 
approximations. Our final result obtained in SD approx- 
imation is 02(5/5/2) = -6.2 ao 3 . 

We calculate the scalar dipole polarizability 00 of Th 3+ 
in 5/5/2 ground state using the expression (Refs. |2(j): 

00(5/5/2) = ^2[Is{nd 3/2 ) + I s (nd 5/2 ) + I s (ng 7/2 )}, (4) 



where 



Is(nlj) 



7 2 

^5/5,2,7113 



9E 



nlj 



E, 



5 /5/2 



(5) 



The breakdown of the contributions to the scalar dipole 
polarizability together with final result for the ground 
state 5/5/2 in Th 3+ are presented in Table 11X1 Again, 
both third-order and all-order results are listed. We use 
the same designations as in Table IxTTlTI We also calcu- 
late the polarizability o corc of the radonlike ionic core in 
Th 3+ . Detailed discussion for the o CO rc in Na, K, Rb, 
Cs, and Fr atomic systems was presented by Safronova 
et al. in Ref. [l^ • We evaluate o corc using random-phase 
approximation (RPA). We find o corc (RPA) to be equal 
to 7.750 (ao) 3 in a.u.. This value was used to obtain 
our final result for the scalar ground state polarizabil- 
ity 00(5/5/2) = 15.1 a 3 ,. We note that unlike the case 
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TABLE VIII: Contributions to tensor polarizability of Th IV in the ground state 5/ 5 / 2 calculated using third-order MBPT 
q (df+2+3) an£ j all-order SD method The third-order and SD dipole matrix elements and corresponding experimental 



transition energies 


are also given. 


All values are 


in a.u. 
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E n — E v 


(DF+2+3) 


7 (SD) 


(DF+2+3) 
"2 


4 SD) 


main 
«2 


5/5/2 


6^3/2 


0.041888 


1.337 


-1.530 


-4.740 


-6.206 




5/5/2 




n C /I £QOQ 

U.040OZO 


n 1 q 
U.zlo 


-u.zoy 


-u.uuy 


-U.U14 




5/5/2 


8(^3 /2 


0.737917 


0.187 


-0.210 


-0.005 


-0.007 




5/5/2 


6d 5/2 


0.066005 


0.362 


0.412 


0.253 


0.326 




5/5/2 


7d 5/2 


0.553263 


-0.036 


-0.052 


0.000 


0.001 




5/5/2 


8^5/2 


0.740531 


0.031 


0.046 


0.000 


0.000 




5/5/2 


557/2 


0.726234 


0.640 


-0.712 


-0.022 


-0.028 




5/5/2 


657/2 


0.827772 


-0.416 


0.413 


-0.008 


-0.008 


<5a!f rc (DF) 












0.042 


0.042 


c4 ail (DF) 












-0.273 


-0.273 


Total 












-4.763 


-6.166 



TABLE IX: Contributions to scalar polarizability of Th IV in the ground state 5/ 5 / 2 calculated using third-order MBPT 
^(df+2+3) a ij. or( j er method ao SD ''- The third-order and SD dipole matrix elements and corresponding experimental 



transition enerj 


;ies are also j 


riven. All values are 


in a.u. 










Contribution 


V 


n 


E n — E v 


7 (DF+2+3) 
^vn 


7 (SD) 


(DF+2+3) 




mam 


5/5/2 
5/5/2 


6d 3 / 2 
6d 5/2 


0.041888 
0.066005 


1.337 
0.362 


-1.530 
0.412 


4.740 
0.221 


6.206 
0.285 




5/5/2 
5/5/2 


557/2 
657/2 


0.726234 
0.827772 


0.640 
-0.416 


-0.712 
0.413 


0.063 
0.023 


0.078 
0.023 




5/5/2 
5/5/2 
5/5/2 
5/5/2 


7^3/2 

8d 3 / 2 
7d 5 / 2 

8^5/2 


0.545323 
0.737917 
0.553263 
0.740531 


0.213 
0.187 
-0.036 
0.031 


-0.259 
-0.210 
-0.052 
0.046 


0.009 
0.005 
0.000 
0.000 


0.014 
0.007 
0.001 
0.000 


a{, ail (DF) 
Q coro (RPA) 
8a c OIC {KPA) 
Total 












0.762 
7.750 
-0.050 
13.523 


0.762 
7.750 
-0.050 
15.073 



of neutral Fr, the core contribution is very large, 50%. 
The calculation of the ground state polarizabilities for 
Fr-likc Th provides another test of the quality of atomic- 
structure calculations. There are no experimental results 
for the Th IV polarizabilities at this time. An accurate 
measurement of the Th IV polarizability, combined with 
these calculations, may be used to derive the values of 
the 5/5/2 — 6d 3 / 2 El matrix elements and to evaluate the 
accuracy of the RPA core value. 



VI. CONCLUSION 

In summary, a systematic relativistic MBPT study of 
energies of ns\/2 {n= 7-10), npj (n = 7-9), ndj (n= 6-8), 
nfj (n= 5-7), and ngj (n= 5-6) states in Fr-like thorium 
is presented. The energies are in good agreement with 
existing experimental energy data and provide a theoret- 
ical reference database for line identification. A system- 



atic all-order SD study of reduced El matrix elements 
and transition rates for the 96 electric-dipole transitions 
in Th 3+ is conducted. Lifetimes are calculated in the SD 
approximation for nlj levels. Multipole matrix elements 
(7si/ 2 -6dj, lsi/2 -5/,-, and 5/5/2 - 5/ 7/2 ) are evalu- 
ated to obtain the lifetimes of the 5/7/2 and 7si/2 states. 
Scalar and tensor polarizabilities for the Th 3+ ground 
state are calculated using relativistic third-order and all- 
order methods. We believe that our energy, lifetime, and 
polarizability results will be useful in analyzing existing 
experimental data and in planning future measurements. 
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